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Abstract 

We study renormalizable Lorentz invariant stable quantum field theories in two space-time 
dimensions with instantaneous causal structure (causal ordering induced by the light "cone" 
time ordering). These models provide a candidate UV-completion of the two-dimensional 

CN ghost condensate. They exhibit a peculiar UV/IR mixing — energies of all excitations become 

(3J[) arbitrarily small at high spatial momenta. We discuss several phenomena associated with this 

mixing. These include the impossibility to reach a thermal equilibrium and metastability of 
all excitations towards decay into short wavelength modes resulting in an indefinite turbulent 
cascade. In spite of the UV/IR mixing in many cases the UV physics can still be decoupled 
from low energy phenomena. However, a patient observer in the Lineland is able to produce 

i— i arbitrarily heavy particles simply by waiting for a long enough time. 

& 1 Introduction 

Our chances for superluminal travel are very slim. There is a single Lorentz invariant causal 
structure in (3 + 1) dimensional Minkowski space-time and it is not compatible with superluminal 
signal propagation. However, life may be more interesting in a (1 + 1) dimensional world — the 
Lineland There one may use a light "cone" coordinate 

(N 
00 

o 



x = t + x 



as an absolute time. Then one may define an unconventional causal structure, such that the causal 
future of any observer at the moment of time x + = a consists of all space-time events in the upper 
right half-plane x + > a, and the causal past consists of all space-time events in the lower left 
^ half-plane x + < a. Lorentz boosts in two dimensions act as 



x + — > Xx + , x — )• A 1 x 

with A > 0, so that the light cone causal ordering is Lorentz invariant. Its existence is related 
to the absence of continuous spatial rotations in the Lineland. All what is left is a discrete 
spatial parity x + -H- x~, which is necessarily broken with this causal ordering (as it is in the 
real world). Following [1] we call this causal structure instantaneous, since it does not admit any 
causally disconnected regions. From now on we will be referring to x + and x~ as time and space 
coordinates correspondingly, unless specified otherwise. 

As shown by [T] , the Einstein-aether theory in two dimensions [21 E] provides an example of a 
renormalizable (in fact, asymptotically free) and apparently consistent theory realizing this causal 
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structural This model is subject to the following subtlety though. At the classical level the 
Einstein-aether theory exhibits a spontaneous breaking of Lorentz symmetry. One may expect 
this symmetry to be restored at the quantum level as prescribed by the Coleman-Mermin- Wagner 
theorem [HE]. However, non-compactness of the Lorentz group complicates direct construction 
of well-defined Lorentz invariant observables in this case. 

Instead, in the current paper we chose to concentrate on simpler examples of instantaneous 
theories such that this subtlety is avoided. Our main goal is to develop a basic intuition about these 
models by performing a number of simple calculations, aimed to demonstrate that instantaneous 
theories are amenable to a straightforward perturbative analysis and possible to work with in spite 
of a number of rather unconventional and, at first sight, counterintuitive properties. 

The origin of these peculiarities is easy to understand. The Lorentz invariant dispersion relation 
in instantaneous theories is 

U = W\' (1) 

where a; is a frequency defined with respect to the time coordinate which is x + and k is a spatial 
momentum with respect to x~ . This dispersion relation exhibits a peculiar form of the UV/IR 
mixing — the higher the spatial momentum of a mode is, the smaller is its energy! In particular, 
the spectrum is always gapless no matter how heavy the "mass" fi 2 of a field is. Understanding 
the consequences of this unconventional dispersion relation will be our main theme. 

The rest of the paper is organized as follows. In section [2] we start by introducing the class 
of theories we will be working with. To keep things as simple as possible we restrict to models 
of a single scalar field <fi. A Lorentz invariant kinetic term with two time derivatives is (d + d-4>) 2 . 
Note that it has four derivatives with respect to the Minkowski time t, so a theory with this 
kinetic term would contain ghosts, if quantized with respect to the conventional causual structure. 
This problem does not arise when quantizing using the light cone time and instantaneous causal 
structure. 

With this kinetic term the naive scaling dimension of the field is —1, so there is an infinite 
number of relevant and marginal operators. For example, any interaction Lagrangian of the form 
P((f>, d + (j)d-(j)) leads to a renormalizable theory according to the naive power counting. In general, 
this will lead to a strongly coupled theory in the IR (i.e., at small uj and k), however, as we will 
see, this strong coupling can be avoided by introducing a sufficiently large "mass" term fj, A (j) 2 . Our 
principal working model of an interacting instantaneous theory will be a 3 theory, which is not 
only super-renormalizable, but also finite. 

To study the physics of instantaneous theories we start by looking at the non-interacting 
model. We calculate the retarded propagator to see explicitly the presence of instantaneous 
signals and discuss an unusual consequence of the dispersion relation ([!]) — the impossibility to 
reach the equilibrium Bose-Einstein distribution at finite temperature. We continue in section 
[3] with introducing the cubic interaction and discuss another interesting consequence of the co/k 
dispersion relation — the instability of all particles towards a two body decay into particles with 
larger spatial momenta. We also argue that in spite of the UV/IR mixing the approximate 
decoupling of heavy particles (i.e., particles with large fi 2 ) is still possible. Nevertheless, the 

1 To be precise, we are talking here only about vector sector of the Einstein-aether, with gravity being decoupled. 
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UV/IR mixing makes high energy experimental physics simpler and less expensive — an arbitrarily 
heavy particle will be eventually produced in the decays of lighter particles, so the discovery of 
the theory of everything is simply a matter of time in the Lineland. 

In order to understand what the absence of the equilibrium Bose-Einstein distribution means, 
in section [4] we use the Boltzmann equation to describe the non-equilibrium behaviour of a homo- 
geneous gas in the instantaneous 3 theory. As a consequence of particle decays an arbitrary initial 
state gives rise to a turbulent cascade towards higher momenta. An arbitrary initial distribution 
function approaches the universal scaling shape at late times. In the limit of small occupation 
numbers we find this scaling solution analytically. Physically, this cascade can be interpreted as 
a process of Bose condensation which cannot be completed in a finite time, because zero energy 
state corresponds to the infinite spatial momentum. We present our conclusions in section [5j 

Before proceeding, a comment is in order. To ensure a better IR (or, to be more precise, 
long distance) behavior it is often convenient to compactify the ^"-coordinate on a circle. One 
may complain that this breaks the Lorentz symmetry. However, from our view-point, the major 
interest of working with a Lorentz invariant theory is that it can be coupled to gravity in a natural 
way (we leave for a future work the study of what happens then). A "soft" breaking of Lorentz 
symmetry, like imposing periodic boundary conditions in x~ , does not spoil this property. 

2 Free Theory 

A simple example of a free instantaneous theory is provided 

The corresponding field equation in the Fourier space is 

(co 2 k 2 - /i 4 )0 = (3) 

so that at any absolute value \k\ of the spatial momentum we find one left-moving and one right- 
moving mode with the UV/IR mixing dispersion relation ([I]). Massless instantaneous theories are 
somewhat degenerate already at the level of a free theory — the dispersion relation ([I]) does not 
give rise to propagating waves for fi = 0. In what follows we restrict to the massive case // 4 > 0. 

In principle, one can also include in (|2| a conventional kinetic term Z /z 2 <9+0<9_0. Then at 
small frequencies and momenta our model would reduce to a conventional scalar field theory 
quantized in the light front. This opens an interesting possibility to interpolate smoothly between 
conventional and instantaneous physics. For simplicity, we restrict to the case Z = 0, this choice 
can be enforced by requiring the invariance under x~ — > —x~ parity. 

A well-known subtlety in the light front quantization is the zero mode problem (see, e.g., [6]). 
The cause of the problem is that if one expands the field (f) into modes with different spatial 
momenta the constant k = mode does not acquire a kinetic term from (pi). 

It turns out this problem ameliorates in instantaneous theories due to the presence of the term 
in ([2]), which is second order in the light cone time. To see this it is convenient to compactify the 
light cone time on a finite interval and to impose zero Dirichlet boundary conditions. This way 



by the following Lagrangian, 



(2) 



3 



the zero mode gets projected from the theory and does not lead to any additional constraints. 
One may recover the non-compact theory by sending the compactification length to infinity 

Note, that this argument cannot be directly applied to conventional theories. First, compact- 
ification on the light-like direction is subtle with a conventional causal structure. Even if one 
ignores this, it is immediate to see that the Dirichlet boundary conditions are incompatible with 
the field equation in the absence of the second derivative kinetic term. Indeed, given some initial 
data at any moment of time one can think of the field equation as an ordinary differential equation 
in x' for the uknown function / = d+(j). This equation is second order for the theory (J2j), 

d 2 J = /iV , 

so that one can impose two Dirichlet boundary conditions. For the conventional theory this 
equation is first order so it is impossible to impose two Dirichlet conditions. 

Clearly, the same reasoning can be applied to project out the zero mode in the interacting 
theory as well. 

Note also that for a conventional theory without higher derivative terms the light cone and 
conventional quantizations give rise to the same results, if done properly. The theory ^ demon- 
strates this is not the case in general. Treated by conventional methods this model contains 
negative energy states — ghosts. As a consequence, already at the classical level, if one adds non- 
linearity into this theory the vacuum becomes unstable, arbitrary small field fluctuations grow 
indefinitely by populating positive and negative energy modes. This does not happen if one uses 
the light cone time to define the Cauchy problem. 

One should not be surprised that different choices of a time coordinate or, equivalently, differ- 
ent ways to pose the Cauchy problem give rise to inequivalent physical results even though the 
Lagrangian is the same. A more familiar example of this kind is a conventional massive theory, 
which turns from being stable into tachyonic if one exchanges the role of space and time. With 
different choices of the Cauchy surfaces one imposes different regularity conditions at the infin- 
ity on the allowed field configuration, and this may tame/induce the instability of the Cauchy 
problem. 

To see explicitly the instantaneous behavior in the theory ^ let us study the field response 
to an external source, which is determined by the retarded propagator, 

G ret (x + ,x-) = — / dkdu—- — = — — ±-J (vVaH>_|) (4) 

47r z J k z u z — fi — leu 2fi z 

We see that there is indeed a non- vanishing field response at all x + > 0. One may worry that 
the fall-off of the retarded propagator at large x~ becomes more and more slow at early times. In 
particular, at x + = + the retarded propagator approaches a non- vanishing constant, indicating 
that the k = infinite frequency mode gets excited. This is an artefact of using the source which 
is switched on and off instantaneously. For physical sources smeared in time (and space) the field 
response goes to zero at the spatial infinity at all times. 

Apart from instantaneous signal propagation the retarded propagator Q exhibits other in- 
teresting features following from the UV/IR mixing dispersion relation ([!]) and persisting in the 
interacting theory. Consider an observer located at a fixed spatial position x~ . After a localized 
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source turns on at the origin, she immediately observes the field response. At late times x + — > oo 
this response asymptotes to 

G re t{x + —> oo,x~) oc (x + \x~ \ ) ~ 1//4 sin(\/ [i 2 x + \x~\) . 

We see that the characteristic wavelength, A ~ /i^ 1 (\x~\/x + ) 1 ^ 2 , of the observed signal be- 
comes shorter at late time. On the other hand the characteristic time-scale for its variation, 
r ~ |) 1//2 , becomes very long. So at late times the observer detects a very short-scale 

and almost static "noise". 

Another piece of information which is usually straightforward to extract from a free theory is 
a finite temperature behavior of a system. However, we find a surprise here. A thermal theory 
has to possess the Bose-Einstein distribution 

< k ) = \ IT > ( 5 ) 

where fl < is a chemical potenital. However, the total particle density given by J n(k)dk diverges 
at large momenta. This is another manifestation of the UV/IR mixing — there is an infinite number 
of high momentum levels with energies arbitrary close to zero, so that one expects to observe a 
never ending process of Bose-Einstein condensation onto the infinite momentum ground state. In 
section [4] we will provide a detailed kinetic description of this process. 

Note, that fermions obeying the UV/IR mixing dispersion relation ([I]) would exhibit the same 
property — changing the sign in the denominator of the distribution ^ does not make its integral 
converging. The physical reason for that is similar — for any finite Fermi momentum there is always 
an infinite number of states below the Fermi surface, so that thermalization cannot be completed. 



3 Interacting Theory 

Let us turn to an interacting theory now. For simplicity, we restrict to the simplest (super)renor- 
malizable interactional 

V int = |0 3 • 


It is immediate to see now the dynamics responsible for the absence of thermalization pointed out 
at the end of the previous Section. Namely, a peculiar property of the dispersion relation ([TJ is 
that two body decays of a particle into a pair of particles of the same mass and larger momenta 
are kinematically allowed. For a particle with a momentum k the momenta of the decay products 
are 

p = —(fk, q = k — p = (f 2 k , 

where if = (\/5 + l)/2 is the golden ratio. Consequently, the momenta of decay products in an 
infinite chain of two body decays originated from a single initial particle can be approximated by 
the Fibonacci numbers. 

2 As usual, the vacuum is unstable in the (/) 3 -theory at the non-perturbative level, but this is irrelevant for our 
perturbative analysis. 
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Let us calculate the corresponding decay rate. The calculation is parallel to what one does in 
conventional theories, one should only be careful to keep track of the correct state and operator 
normalizations. We introduce canonical creation and anihilation operators a k ,ak satisfying 

[a k ,al,] = 6(k - k') . 

For and and its canonically conjugate momentum tt^ to satisfy the canonical commutation 
relations, the field operator should be defined as follows 

w* „vw\ ab+cx ' (6) 

Note the extra factor in the wave function normalization compared to the standard (2w) _1 / 2 
[7J, it is related to the presence of two spatial derivative in the expression for the canonical 
momentum, tt^ = — <9 + <9^0. From here the standard expressions for the Feynmann rules and 
transition rates [7J adopted for the theory at hand give the following result for the two-body decay 
width, 

r » 2 A, , W-p-gKfi-gr-fa) f ( iV 1 m 
r = 8Si irt N = w(^%J (7) 

Note that the ratio of the width to the frequency of a particle V/u does not depend on the spatial 
momentum k. This behavior could have been anticipated from the Lorentz invariance of the 
theory, so that the result ^ for the width T is fixed by the Lorentz invariance and dimensional 
analysis up to an overall numerical factor. As a consequence, even though the width itself diverges 
at small spatial momenta, particles remain narrow resonances and the theory is weakly coupled 
(at least, as far as this particular process goes), provided the dimensionless coupling gf fj, 5 is small. 

It is natural to expect that the dispersion relation ([T| makes decoupling of heavy particles 
subtle compared to the standard case. As an illustration consider a situation when in addition 
to a light field of mass /i, there is another one with mass M ^> p. Normally, at energies below 
M heavy particles simply never get produced. Consequently, one can write an effective low- 
energy theory of light fields only, where all the memory about heavy particles present in the full 
microscopic theory is absorbed into the values of coupling constants describing local interactions 
of the light fields. This is no longer the case with the UV/IR mixing dispersion relation ([!]). For 
example, a two body decay 

A;) -> \ti,p)\M,q) (8) 

with \k\ ~ fi and \p\, \q\ ~ kM 2 /fi 2 is kinematically allowed no matter how heavy the mass M is. 
Consequently, in a strict sense, the full decoupling never takes place and one can always produce 
arbitrarily heavy particle simply by waiting long enough. 

However, if the coupling constants for light and heavy fields are of the same order, there is an 
approximate notion of decoupling. One indication for this comes from the possibility to perform 
a Wick rotation, which for instantaneous theories amounts to replacing x + — > ix + . The resulting 
"Euclidean" propagator (co 2 k 2 + /i 4 ) 1 does not have poles just like it happens in ordinary theories. 
Similarly to the conventional case for calculating Green's functions involving light fields only at 
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low values of kinematic invariants one can Taylor expand the heavy propagators in o; 2 /c 2 /M 4 , 
which trades a heavy field for an infinite number of local vertices involving light fields only — the 
standard "integrating out" procedure. 

For low energy processes involving on-shell heavy particles, such as the two-body decay ([§]), 
the suppression comes from the difference in the wave function normalization in ^ between heavy 
and light field^j Replacing an external line of a light field by a heavy one in a given Feynmann 
diagram brings in a suppression both because of the heavier mass in the denominator of ([6]), and 
because the corresponding spatial momentum is higher as the example (|8j) illustrates. 

To further illustrate this reasoning let us consider a super-renormalizable (and finite) theory 
of three scalar fields <fii, (f>2 and 4>m with masses fi, /i, and M satisfying the hierarchy fi <C M. As 
an interaction potential we chose 

Vim = 912(f) M (01 + 02) + g((pl + <j>l) , 

so that different light fields interact only through the exchange of the heavy field. Note first that 
the decay of a light particle into two heavy particles has a width 

J- Ihh 



fi 2 M e k ' 

where we suppressed all order one factors. Consequently, in order for this process to be subleading 
with respect to the processes involving light fields only one needs the condition 

9i2 < g—r 9 

to be satisfied. 

In particular, as anticipated, processes with heavy particles are indeed suppressed if gi2 ~ g. 
In this case, just like in conventional theories, we can integrate out the heavy field, which will 
result in direct interactions between different light fields 0i, 02- At the leading order in g 12 these 
are mass, kinetic and higher order in momenta mixings of the form 

gl.M-^Md+d-) 2 ^ , (10) 

which are small as a consequence of ^ and the smallness of g//J? . 

To summarize, we see that in many situations decoupling in instantaneous theories works 
similarly to the standard case, with the important difference that with enough time one will 
always be able to produce arbitrarily heavy particles on-shell. On the other hand, the condition 
^ is not necessary for the theory to be perturbative (for instance, one can set g — 0). If it is 
violated, one ends up being in an interesting situation when the dominant processes involving 
low-energy light particles are decays into heavy particles. 



3 This suppression is not unrelated to the Euclidean argument above. At the end of the day, the propagator is 
built as a convolution of two wave functions. 
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4 Turbulence 

After developing some experience with simplest interacting instantaneous theories, let us come 
back to the observation at the end of Section [2] that thermalization is impossible in these models 
and describe quantitatively the turbulent cascade corresponding to the Bose condensation onto 
an infinite momentum ground state, which was anticipated there. 

The first simple thing to try is just to study numerically the evolution of some smooth ini- 
tial data in the interacting 3 theory The typical result is shown in Fig. [!} As expected, we 
indeed observe that at late time the field configuration turns into a short-distance noise. To de- 
scribe the statistical properties of this noise quantitatively we make use of the kinetic approach. 
Namely, we consider a homogeneous (on average) gas of 0-particles with spatial momenta distri- 
bution described by the distribution function n(k). The evolution of this distribution function is 
determined by the Boltzmann equation, where in the limit of weak non-linearities the collision 
integral contains two-body decays, and the inverse process — two-to-one collisions. A particle with 
a momentum k participates in six processes 

k -H- —{pk,ip 2 k ip~ 2 k -H- k, —ip~ l k — ip~ l k -H- k, —ipk (11) 

The k — )■ pq decay contributes to the collision integral as 

'•' f <p + -) n(k)(n(p) + l)(n(q) + l) (12) 



and the inverse processes as 



8fi 8 k \ (p 



<J ' ip+ - ) {n(k) + l)n(p)n(q) . (13) 



fj, 8 k \ (p 
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The resulting Boltzmann equation reads 

d T n(k,r) = — f — n(k,r) + n(— <^ _1 A;, r) + n(ip~ 2 k,r) — 
k\ 

—n{ip 2 k, r)n{k, r) + n(ip 2 k, r)n(— (p~ 1 k, r) — 2n(—<fk, r)n(k, r) + 

+n(Lf~ 2 k, r)n(k, r) + n(ip~ 2 k, T)n(—Lp~~ 1 k, r) + n(—Lf~ 1 k, r)n(—(pk, r)J (14) 

where we absorbed the dependence on the coupling constant into the definition of the time variable 

i g 2 ( . IV 1 



27r8/i 8 V <P 

Let us start by considering a low density gas with all occupation numbers being smalQ n(k) <C 
1. In this regime one can drop quadratic contributions to the collision term in ( Jl4| and all the 
dynamics is determined just by the particle decays, 

d T n(k, t) = \ (— n(k, r) + n(— tp~ l k, r) + n(ip~ 2 k, r)) . (15) 
k 

This equation turns out to be amenable to the analytic study. For simplicity we restrict to the 
case n(k) = n(—k). It is natural to look for a special solution of (15) of the form 

n (k,r) = f (log-) . (16) 



Then the function f(x) satisfies the following equation 

- e x f'(x) = -f(x) + f{x - \ogip) + f(x - 2 log <f) . (17) 

By performing the Fourier transform we reduce the problem to finding an analytic function de- 
caying along the real axis and satisfying the following quasi-periodicity condition 

{-iy - l)f(y - i) = (1 - ^ +1 )(1 + ^)f{y) ■ (18) 



In the absence of the cp dependent prefactors on the r.h.s. of (18) this would be a defining property 
of the T-function, T(—iy). The extra prefactors are straightforward to account for by writing the 
solution in the form 

f(y) = T(-iy)P + (y)P-(y), 
where P± are two uniformly converging infinite products, 

oo oo 

p + ( y ) = n (i - ^ y - m ) , p-{ y ) = n t 1 + v w - m ) (19) 

m=0 m=2 

which can also be expressed through the q-Pochhammer symbols as P± = (±(p ±1+iy , <^) 00 /(l =p 
tp ±1+iy ). Note, that all the poles of the Gamma function get cancelled by the zeroes of P + . Gamma 



1 Note, that this is a quantum regime, which cannot be captured by solving the classical field equations. 
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Figure 2: Time evolution of the distribution function as a function of log A; for the linearized (left) 
and full (right) Boltzmann equation. In both cases the distribution function converges to the 
universal shape at late times, which in the linearized case well agrees with the analytic solution. 



function also guarantees the convergense of the inverse Fourier transform at the infinity. This all 
implies that the function 

f(x) = J e-^T(-iy)P + (y)P^y)dy (20) 



satisfies the equation (17) and decays at the infinity faster than e a ' x ' for any a. Clearly this 



function is real, however in order to correspond to a physical solution of (17) it also has to be 



positive. Fortunately, as the plot in Fig. ^a) shows, f(x) is indeed positive and, consequently, 



provides a solution to the linearized Boltzmann equation via (16). Moreover, the same plot 
demonstrates that this solution is actually a universal attractor, describing the shape of the 
distribution function as a function of log k at late stages of the cascade independently of the 
initial shape (of course, there is a freedom in choosing an overall normalization, which is fixed by 
the total energy of the gas). 

To understand the origin of this attractor behavior note that by applying constant time shifts 
to uq we can construct the whole family of solutions of the form n u (k, r) = f(logk — log(r + u)). 
All these solutions have the same bell-like shape as a function of log k, and differ only by an 
overall shift. Given an arbitrary initial profile n^k) at r = we can then decompose it in a linear 
combination of n u (k, 0), with u > 0. The corresponding solution of the Boltzmann equation takes 
then the following form 

POO 

n(k,r)= / duf (log k - log (r + u))c(u) (21) 
Jo 

where c(u) are the coefficients in the linear decomposition of rii(k). For the initial configuration 
to have a finite total density of particles the function c(u) should decay at the infinity faster than 
1/u. Then the late time asymptotics of the solution is 

n(k, r — )• oo) = / ^log —J J duc{u) 
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which is exactly the attractor behavior exhibited in Fig. |2 



The same argument cannot be applied for the full non-linear Boltzmann equation (14). How- 



ever, as illustrated in Fig. |2p), numerical results indicate that the attractor behavior persists at the 
non-linear level as well. Independently of the initial shape at late time the distribution function 
is described by a universal scaling solution f n i(\og -,E), where E is a total energy. Heuristically, 
the origin of this attractor behavior can be understood as a consequence of the rate of interactions 
getting lower as the momenta of the particles grow. As a consequence, for any initial distribution, 
"fast" low momenta particles always have enough time to catch up with "slow" high momenta 
modes and to establish a stable shape, which keeps drifting into UV afterwards. 



5 Discussion and future directions 

To conclude, we hope to have convinced the reader that instantaneous theories present a tractable 
example of apparently local Lorentz invariant quantum field theories exhibiting interesting non- 
local features. Our analysis of their dynamics is far from being exhaustive, we focussed on the 
simplest case of two-body decays. It would be interesting to study more general processes, in 
particular one may expect peculiar IR singularities in scattering processes following from the 
absence of the gap in the spectra of instantaneous fields. 

As we said in the Introduction, the main motivation to insist on Lorentz invariance is that it 
provides a natural way to couple these models to gravity. It will be interesting to study whether 
consistent instantaneous quantum gravity models can be constructed and what is the fate of black 
holes in these theories. 

Note also, that a two dimensional version of the ghost condensate model [8] becomes power 
counting renormalizable if promoted to an instantaneous theory 

£ 9C = (D0) 2 + P((«90) 2 ) . (22) 

It will be interesting to see whether this setup provides a useful lesson for understanding the 
dynamics of the ghost condensate and more general phases of massive gravity [9]. Dynamics of 
these theories may be somewhat more subtle than in examples studied in this paper, because the 
shift symmetry acting on prevents one from adding a mass term /i 4 2 . 

As a totally different direction, note that one may view the theories studied here as a specific 
case of non-relativistic theories with an anisotropic scaling symmetry 

r — > At, cr — > \ z cr. 

Theories of this kind with z < are likely to share many of the features of z = — 1 case studied 
here. It is an interesting question whether anisotropic scaling symmetries with z < may arise 
in condensed matter system. In this regard it is worth noting that the dispersion relation 
describes short-wavelength nondivergent Rossby waves in the Earth atmosphere in the absence of 
mean currents [10J. 

Finally, the physics of instantaneous models relies crucially on special properties of a two- 
dimensional Lorentz group. It is not clear whether the lessons learnt here may be applied in 
higher dimensions. A promising proposal in this direction [TJ may be to study the dynamics 
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of a string, with the instantaneous causal structure on its world-sheet, propagating in a higher 
dimensional space-time. We hope to report a progress on this and other related topics in a near 
future. 
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